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Abstract: 

Whether our universe with present day acceleration can eventually collapse is 
very interesting problem. We are also interesting in such problems, whether the 
universe is closed? Why it is so flat? How long to expend a period for a cycle of 
the universe? In this paper a simple "slow-fast" type of the quintessence potential 
is designed for the closed universe to realize the present acceleration of our universe 
and its eventual big collapse. A detail numerical simulation of the universe evolution 
demonstrates that it divides the seven stages, a very rich story. It is unexpected that 
the quintessential kinetic energy is dominated in the shrink stage of the universe with 
very rapid velocity of the energy increasing. A complete analytic analysis is given 
for each stage. Some very interesting new problems brought by this collapse are 
discussed. Therefore our model avoids naturally the future event horizon problem 
of the present accelerative expanding universe and maybe realize the infinite cycles 
of the universe, which supplies a mechanism to use naturally the anthropic principle. 
This paper shows that the understanding on the essence of the cosmological constant 
should contain a richer content. 



I. INTRODUCTION 



The birth and termination of the universe is a riddle for mankind. An old story 
is that a matter-dominated closed universe will shrink finally t^' , so that the universe 
begins its new big bang through a bounce, and repetition forever to realize infinite 
kinds of universes. However, this old idea was challenged by a fact known recently, 
i.e., our present universe is accelerated expandingPl. If the reason to cause the 
acceleration of the universe is the real cosmological constant, then the universe will 
be accelerated expanding forever, never shrink in its future, thus we can not see the 
cycle of the universe. It is believed generally that the appearance of the wisdom life 
needs some kind of adjustment of some physics parameters, so that we need many 
worlds to realize this adjustment, that is the anthropic principle '^1. If the universe 
never shrink, we have to invoke the concept of the chaotic inflation'^' . The universe 
is inflating eternally in a larger region, produces many many different worlds in a 
fractal manner. Today this new idea is very popular in the cosmologist community. 
Comparing earnestly the old and new ideas, one can find out that the argument 
of the chaotic inflationary production of many universes owns some difficulty of 
mathematical description, and the old idea is more simple and clear. 

We hope that our universe is accelerative in present and will shrink finally. The 
advantage of this idea is that the universe maybe stays in its infinite loop to realize 
many world scenario, that this model agrees with the present decisive astronomi- 
cal observations, specially the cosmic microwave background radiation (CMBR)'^', 
and that it avoid a difficulty about an event horizon for a perpetual accelerated 
universe'^l. In order for the universe to shrink finally there must be some negative 
energy in it, for instance, a negative cosmological constant or a positive-curvature 
negative-energy density. In the first case we need that the cosmological constant 
changes from the present positive value to the future negative one, thus the cos- 
mological constant is not a true "constant" at all. It is better for us to accept the 
quintessence scheme'^', which has a good basis of the field theory. Anyhow we hope 
to avoid a "true cosmological constant" further, no matter what it is positive or 
negative, which is a big obstacle to construct a string (M-) theory'^l Therefore we 
take the vacuum of the quintessence, i.e., the minimum of its potential, to be the 
zero point of the energy, that is the true cosmological constant is zero absolutely for 
some unknown reason'^'. The force to shrink the universe is remained only as the 
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positive-curvature negative-energy density, i.e., the universe is closed. In spite of 
the universe is closed, we shall show in the section 10 that a reasonable assumption 
from inflation of the early universe gives a very small non-flatness of the universe, 
which is unobservable at present. 

It is a very interesting problem how to shrink finally for a closed universe, which 
has acceleration at present caused by a special quintessence. It should be known 
that for a generic quintessence, for instance with a pure inverse power law potential 
or a pure exponential one, it can not arrive the above goal. That the universe is 
accelerative at present and will shrink in future is to require that quintessence is a 
"(early) slow - (late) fast (rolling)" type. However, the inverse power quintessence 
is a "fast-slow" type, and the exponential quintessence is a "slow-slow" or "fast- 
fast" ones. So that we must design a reasonable "slow-fast" quintessence. In this 
paper we do it in the section 2. This potential must have as simpler form and fewer 
parameters as possible. This special potential is called by us as "Niagara" one, since 
its configuration is very flat for its negative fleld value and falls very rapidly to the 
zero for its positive fleld value, and the quintessence rolls down from a high potential 
hill to a low valley without boundary. 

Whether there is some new story about the evolution of the universe with our 
"Niagara potential"? Yes, some rich and wonderful scenario appears unexpectedly. 
As a typical process from the Big Bang Neucleonsythesis (BBN) to the big collapse 
of the universe, the evolution may be divided into seven stages, as described in the 
section 3. The flrst and second stages are the radiation and matter dominated epochs 
respectively, which are familiar for us. The third stage is the quintessence potential 
energy dominated, and the universe stays in its accelerative phase, in which the 
quintessence rolls slowly with its negative fleld value (see analysis in the sections 4-5). 
However, the quintessence rolls inevitably to its positive field value. When it passes 
through its turning point, the quintessence potential energy decreases suddenly, thus 
the quintessence kinetic energy will be dominated, and this is the fourth stage of 
the universe evolution, which is decelerated (see analysis in the section 6). Since 
the evolution law of the quintessence kinetic energy is faster than matter, when the 
cosmic scale factor expand to some value the quintessence kinetic energy will be 
approximately equal to the matter energy, the universe comes into its fifth phase, 
i.e., the "matter leading attractor" one, since its evolution law is like of matter and 
this new term will be explained in the section 7. Moreover, the evolution law of 



the matter energy is faster than the curvature energy, at some time the universe 
is dominated by the curvature energy, i.e., the sixth phase of the evolution (see 
analysis in the section 8). When the curvature negative energy cancels just all of 
other energies, the universe stop its expansion, and begin to shrink, i.e., the seventh 
phase. We note interestedly that when the universe begin to shrink, the quintessence 
kinetic energy, the matter one and the curvature one are all approximately equal. 
As the universe shrinks, the quintessence kinetic energy increase more quickly than 
other form energies, soon after it the quintessence kinetic energy is dominated in 
the collapse universe, therefore the seventh phase of the universe evolution is the 
"quintessence kinetic energy collapse", it is very different from the old concept of 
the matter collapse (see analysis in the section 9). This will bring many interesting 
problems to deserve further investigation (see the section 11). 

It is worth to be emphasized that it is very necessary to finish a numerical 
simulation in computer for a whole evolution process of the universe, otherwise 
we can not supply a time-scale concept for each stage and a diminishable scenario 
of the universe. We noted often that many seeming feasible models are actually 
impracticable, specially to obtain a suitable accclcrativc period and a reasonable 
age for the present universe. We give a detail result of the simulation of this seven 
stage evolution in this paper. We need not only the detail simulation (in the section 
3), but also an analytic analysis for its essence in a manner of the quantitative 
estimation (in the sections 4-9). 

The model parameters taken by us is only a special point of the wide parameter 
space in our model, however it is typical for an understanding of the essence of the 
evolution. To explore the various behavior of the whole parameter space is not the 
task of this paper. We only hope to express an idea about the present accelerated 
expanding universe still has a chance to shrink eventually. This opens a way to 
realize infinitely many universes. This kind of many universe is more simple and 
easy to be described in order to serve the anthropic principle. Finally we make a 
simple conclusion in the section 12, specially about interesting problems worth to 
be investigated further. 
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II. THE POTENTIAL OF THE QUINTESSENCE 



We have announced that it is necessary to replace a true cosmological constant 
with the quintessence field 4>. The quintessence potential we shall take in our model 
has the following form 

'^f^' - l + expW/M,) - W 

When (f) — > — oo, the potential V{(f)) approaches to Mq , which emulates reahstically 
a true cosmological constant, we take Mq = 0.95 x 10~^°Mp, here Mp = (SttG)"^/^ 
= 2.4 X lO^^GeV is the Planck energy scale. In other hand when (f) is larger than 
zero, the potential V{(f)) declines exponentially and approaches to zero quickly. The 
slow roll parameter of this potential is 

2 (1 + e-^^/^p)2' ^ ^ 

The parameter A will control the velocity of roUing. This potential is "slow-slow" 
type for A < y/2 and "slow-fast" one for A > -\/2. In our model we take A = 2 as 
a typical analysis. Due to the shape of this potential likes the Niagara fall, we call 
it as "Niagara potential" . The field value = 0.44Mp is the turning point between 
the slow and fast rolling. When < 0.44Mp the quintessence is slow roll with e < 1 
and when > 0A4Mp it rolls fast with e > 1. The model with A = 2 is a typical 
"slow-fast" potential. 

The potential of the inverse power law Mo"'""0~° of the quintessence is the "fast- 
slow" type which can realize the interesting tracker solution^^^l . The pure exponential 
potential exp(A0/Mp) is "slow-slow" one for A < and "fast-fast" for A > v^. 
Of course we can design a more compound "fast-slow-fast" potential to realize both 
tracker and shrinkage, but this is not a purpose of this paper. 

III. THE SIMULATION OF THE UNIVERSE EVOLUTION 

The whole evolution obeys the following famous equation, which was researched 
for thousands times in the literatures, 

= ^^{P<P + Pm + Pr+ Pc), (3) 



+ 3H<j) + y = 0, 



(4) 



where H = a/a is the Hubble parameter and a is the cosmic scale factor, t is the 
cosmic time, Ptj, = Pk + Pv is the quintessence (total) energy density, in it pk= 
0^/2 is the quintessence kinetic one, p^, =1^(0) is the quintessence potential one, 
Pm = Pmoici/do)'^ is the matter one, p^ = ProicL/do)'^ is the radiation one, and 
Pc = Pcq{0'/ O'q) '^ is the curvature one which is negative for a closed universe. We 
noted pmo, Pro, Pco are the parameters of the model together with A and Mq, total 
five. A set of initial condition, ao, 0o, 0o at the initial time to, determines uniquely 
the whole evolution from the BBN to final constriction of the universe. As a typical 
analysis, we take the following parameters, to = lsec= 0.4 x 10^^ for BBN time, 
PrQ = 2.6 X 10~^^Mp, pmo = 2.2 X 10~^^Mp for an almost realistic universe model. 
We shall see that it gives a set of the suitable parameters of the present universe. 
The initial condition is taken as 00 = —2 in order to have an enough slow rolling, 
and 00 = 0, i.e., no roll at beginning, we shall see that this assumption is reasonable. 
As for the curvature energy, we take pco = — lO^^^^Mp, we shall explain why this 
value should have such a small magnitude orders, ao is arbitrary unit for a while, we 
shall determine which value it should have. If we know the values of a, and (it 
is determined by function 0(t)) at an arbitrary time t, we can calculate the values of 
all other 15 quantities at this time, such as H, p^, Pr, Pc, Pk, Pv, the various density 
ratios Qj = Pj/pcrit, where the critical density is pcrit = SM^H^, the quintessence 
state equation = (p^ — Pv)/{pk + Pv) and the acceleration parameter q = ad /a? 
in according to the formulae just given, as well as 

a= ^^{pv-'^Pk- Pm/'^- Pr)- (5) 

The key is to obtain a, and as the functions of time. Remember that this 
principle can guide our later analytic analysis. 

It is impossible to obtain the exact solution of the whole process of this set of 
equations Eqs.(0) and (^. We do a fine numerical simulation. The result is unex- 
pectedly rich. The evolution divides seven stages, all listed numbers are obtained by 
computer calculation, the purpose that some numbers preserve their high precision 
form is for convenience in order to compare with the later analytic analysis. 

The result of the simulation is mainly given in the Figs. 1-4, the evolutions of g, 
a, and 0. Some quantities have too large magnitude orders, so that we have to 
use a double logarithmic coordinates. The seven stages of the evolution can clearly 
be seen in the figure of q value, as a feature quantity of the different stage. It 



is a pity that wc arc not able to put the wonderful evolution figures of all other 
14 quantities in this paper for lack of a page space. Due to the continuity of the 
evolution process, the turning points among the different stages have to be chosen 
carefully, even subtly, otherwise it is not easy to obtain the good results in order to 
compare with the analytic analysis. 

1) The radiation dominated epoch, from to to ti — 0.6 x 10^^ M~^, where time 
ti is defined as pri = Pmi, i-e., the radiation- matter equal. The cosmic scale fac- 
tor increases to ai = 1.2 x IO^Oq. The quintessence field value is almost invariant 
(pi — —2Mp. The quintessence velocity is 0i = lO^'^^'^M^. The acceleration param- 
eter takes its value g(o^i) — —1 and qi = —0.75, i.e., the universe is decelerated 
expanding in this stage. The quintessence state equation is almost w^q^i) ~ — 1. 

2) The matter dominated epoch, from ti to t2 — 0.86 x 10^° M~^, where time 
t2 is defined as = '^Pv2, i-e., the matter- "cosmological constant" approximately 
equal. The reason why we choose a coefficient 2 is that at this time the universe 
comes just into the accelerative phase, q2 — 0.005. Before this point q{i^2) — —0.5. 
The quintessence has 02 = -1.997Mp, 02 = lO-^^-^M^, w^2 = -0.9999. The cosmic 
scale factor increases to 02 = 2.4 x lO^ao- 

3) The "cosmological constant" dominated epoch, from t2 to ^ = 3.2 x 10^^ M~^, 
where time is defined as 0(^3) —0, i.e., the ending of the quintessence slow rolling. 
In this epoch the universe stays in the infiationary phase ^(2-^3) — 1 and qs — 0.676. 
The quintessence has 03 = 0.026Mp, 02 = lO'^^'^M^, w^s = -0.784 and ty<^(2-»3) ^ 
— 1. The cosmic scale factor increases to 03 = 1.7 x IO^^Oq. Here the result of the 
point 3 is new. 

In particular, this epoch contains a special point, our present universe, denoted 
as 'V. We have tu = 1.56 x 10^° = 13Gyr, i.e., the universe age. Maybe you 
noted that this number solves the crisis of the universe age'^^l due to the existence 
of the pseudo "cosmological constant". The time point t„ is defined as Pt,(t„) = 
(7/3)pm(^u); i-e., ^mu = 0.3 and flyu = 0.7. The present universe entered in the 
accelerative phase not long ago, g„ = 0.55. In this point we have pru = 10~^^^Mp 
which is just the energy density of the present CMBR with — 2.726°i^. The 
Hubble constant is = lOO/ikm-sec-Mpc"^ = 0.62 x 10"^° Mp , then h = 0.7. 
Through this data we can see the model parameters are chosen reasonably. The 
scale factor of the present universe is a„ = 4 x IO^Oq. 

4) The following results are new, you must note the jumping changes of various 



parameters. The fourth stage is the "kinetic energy leading attractor" epoch, from ^3 
to t4 — 1.6x 10^°^M~^, where time ^4 is defined as PmiU) —0.3pv{t4), i.e., the matter- 
quintessence approximated equal, why such choice will be explained in later analysis. 
After ^3 the universe turns suddenly into the decelerated phase, ^(3-^4) ~ — 1 and 
54 = -0.84. The quintessence has 04 = 95AMp, 04 = lO^^^^-^M^, = 0.25 and 
^<l>{3^4) — 1/3. The cosmic scale factor increase to 04 = 1.8 x IO^^Oq. 

5) The "matter leading attractor" epoch, from ^4 to U = 0.28 x lO^^^M^-^ 
where time is defined as (^^(ts) ~ —1, i.e., the curvature dominated. The scale 
factor of the universe is 05 = 3.1 x During this epoch the universe is in the 
decelerated phase, ^(4^5) — —0.5 and ^5 = —1.55. The quintessence has 05 = 280Mp, 

05 = 10~^®^-^Mp, — 0.26 and iy</,(4^5) ~ 0. The density ratios of the simulation 
are flkb = 0.823, ftrnb = 0.770, Q„5 = 0.482 and Qcb = -1-075. 

6) The "curvature dominated" epoch, from ^5 to — 0.85 x lO^^^M"^, where 
time ^6 is defined as H{tQ) —0, i.e., the absolute value of the curvature energy 
equals the sum of all other energies, called it as the "balanced point" , and the scale 
factor of the universe arrive its maximum value ag = 4.5 x IO^^Oq. During this 
epoch the universe is still in the decelerated phase Qq —00. The quintessence has 

06 = 281. 5Mp, 06 = lO-^^^-^'^M^, w^Q = 0.805. In fact the computer value of Hq 
is not true zero due to lack of the calculation precision. In this point we have the 
interesting proportions Pke/lPcel = 0.461, Pme/lPcel = 0.489 and p^e/IPcel = 0.050 
for various energies. We avoid to use here due to their infinity at the balance 
point. 

7) The big shrink epoch, from to ~ 1.74^6 for a rough estimation, where 
time ts is defined as Og — 0, i.e., the cosmic scale factor shrinks suddenly to almost 
zero! However it is very difficult to look for this point tg in computer, instead of 
we choose another neighboring point tj = 1.47 x lO^^^M^"^ = 1.73^6 to study the 
behavior of the cosmic constriction. At this point ay = 0.27 x IO^^Oq — 0-6/17, 
indeed the universe is shrinking. We have 07 = 289. IMp, 07 = 10~^^^Mp ~ 470006, 
w^7 ~ 1, i.e., the quintessence kinetic energy becomes large rapidly. The computer 
calculation gives 0^7 = 0.99981 and iljnr — 0.00021, 0,^7 — —0.00003, the rests are 
almost zero, i.e., it is a "quintessence kinetic energy dominated shrink". This is a 
very interesting new result obtained by this paper! We shall illustrate its important 
meaning in later. 

Maybe one notes that many terms and scenario have not been illustrated in 
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above text. All things can be done in the following detail analysis, we need patience. 
This numerical simulation demonstrates a vivid story of the universe evolution. In 
order to make one believe this simulation, we must do some convincing analytic 
analysis, and look for the profound relations between various data. We must answer 
the question such as what is the features of the balanced point, the essence and 
inevitability of the shrink, and the reason of the closed universe with sufficient 
flatness. 



IV. THE EVOLUTION OF THE QUINTESSENCE IN THE RADIATION 

OR MATTER DOMINATED 

At the first and second epochs we meet a problem how the quintessence evolves in 
the radiation or matter dominated environment, pn ^ p<^ ^ |pc|, where pn = PnoO'~"' 
with n = 3 for matter and n = 4 for radiation. Suppose the quintessence is in slow 
roll process, i.e., < and 0^0. We have a oc t^/" and H = 2/{nt), and from 
Eq.®, 

3.l0_^ L^. = 0. (6) 

Integrating it we obtain the time te at which the quintessence ends its slow rolling 
from the initial value 0o to the final one 0e — (here index "e" means "end"). 

It doesn't matter for this estimation to use 0e — rather 0e — 0.44Mp. The time te 
must be larger than the present universe age t^. Otherwise, there is not an enough 
time to allow a phase in which the cosmological constant be dominated, since the 
quintessence already rolls fast. Once the quintessence comes into fast rolling, the 
evolutionary law of the quintessence will be changed, does not obey this formula 
rather the later Eq.(|l2). From it we know that it is safe to take 0o = ^2, since 
te = ^Atu- From Eq.(^ we also get the quintessence evolution of slow rolling in 
the matter or radiation dominated epochs 

~ -X-^Mp\n{e-''^''^'^ - nX^M-^M^ty 12) = (f)_{t), (8) 

moreover which is also valid in the whole epoch of the quintessential slow rolling, 
even if the quintessence dominated, i.e., to < t < t-^. Using Eq.(||) we can estimate 



01 = -2(1 - O{10-^^))Mp, 02 = -1.997Mp and 0i~02^ lO'^^M^, coincides with 
the simulation value. Hereafter the wave "~" means the analytic estimated values, 
the estimation of the theoretical values of ai and a2 is elementary and familiar, and 

the values of other various quantities at the time point 2 are all get in according to 
the mentioned principle, thus we omitted these calculations. 



V. THE SLOW ROLLING QUINTESSENCE DOMINATED 

At the third epoch we meet a problem how the quintessence evolves in the 
quintessence dominated, ^ p„ ^ \pc\, and e < 1. This problem has been 
solved partly in the last section, here we concern the scale factor. This is a typical 
inflationary process a oc e^* with H ~const., in which the quintessence rolls slowly 
from (pb — 02 to 0e < OAAMp (The index "6" means "begin" in this paper), we can 
use its standard formulas, i.e., the inflationary e- folding N is given by 

a r't'^ n V 
N^ln-^ M;2(--)#. (9) 

Using it we have 

N = A-^M-i(0e - ^e-^^«/^-) + A-iM-i(-06 + ^g"^^"/^-), (10) 
A A 

and the time at which the quintessence slow rolling is terminated, 

te(3) ~ V3X~^MpMQ^e'^'^>'^^^, (11) 

which must be larger than the present universe age t^- 

However we are not able to analysis the data of the time point 3 for a while since 

it involves the next stage evolution. After we finish the analytic analysis of the next 
stage, we shall return back to do a theoretical estimation of the data of the third 
cosmic stage. 



VI. THE QUINTESSENCE KINETIC ATTRACTOR 

A question we face in the fourth stage of the universe evolution is the fast rolling 
quintessence dominated, i.e., from the time points 3 to 4, ^ p„ ^ \pc\, with 
> Mp and e > 1. It is unexpected that its solution is very simple 



= 2X~^Mp \n{XM-^M^t) = (12) 

which can be easily checked by putting it into Eqs.(^) and (|), furthermore only this 
formula is still valid in the fifth and sixth stages, other formulae will different. Its 
character is that the quintessence state equation is = X^/S—l and the evolution of 
the cosmic scale factor is a oc t^-^ ^ ^ the acceleration parameter g = 1 — X'^/2. When 
we take A = 2 the evolution of this quintessence is just like radiation Wr = 1/3, but 
it is only a fortuity. In this case its kinetic energy is twice of its potential, pk = 2p^. 
It is seen clearly that there is not the decelerating of the universe for the "Niagara" 
quintessence with A < v^. 

Now a careful observation is needed for the slow-fast turning point. When the 
quintessence rolls slowly it obeys the evolution Eq.(^, and when it rolls fast it obeys 
the evolution Eq.(|l2]). Observing the neighbor of the point 3 of the simulation curve 
in the Fig. 3, which has a subtle transition, we have a good approximation about 
this point. We suppose that = 0_ for t < and = 0+ for t > and 0- = 0+ 
for t = t^. These two curves intersect at the time point ts, we can obtain the 
cross point k = 2.7 x lO^^M'^ and 0_( tg) = 0+( h) = 3.888Mp by solving union 
equations. This point is near the slow-fast turning point, = 0.44Mp, therefore the 
inflationary e- folding formula Eq.(|lOD is still valid, we take the slow rolling ending 
point as 0e =0-( ^s) and the beginning point is (pb — — 2Mp then we have from 
Eq. ([lO|) = 16.5, then = 1.5 x 10^52- Due to the delicacy of the transition point 
3, it is reasonable to take the theoretical value of 03 as almost zero (therefore different 
from the value 0_( ta), a skill for us), and obtain 03~ 10~^^ Mp from Eq. (|l2D . These 
estimations are good coincident with the simulation data. 



VII. QUINTESSENCE-MATTER ATTRACTOR 

A question we face in the fifth stage of the universe evolution is the fast rolling 
quintessence domination with the matter leading, i.e., from the time points 4 to 5, 

^ Pr ^ \pc\. Putting a (X t^/" and H = 2/{nt) into Eqs.(||) and (^), we 

can find out that 

n 1 ?) n 3 

n^ = i-- = - n, = — = - = ___ = _. (13) 

This the attractor solution is given by Ref.[12], we call it as the "matter leading at- 
tractor" if n = 3 and the "radiation leading" one for n = 4 due to their evolution be- 



havior. Now we can understand why we choose the time point 4 as Pmii^) =0.3pt,(t4) 
since the matter density increases from Vtm/^v — to Vtm/^v — 2/3 when the uni- 
verse comes from the quintessence leading phase to the matter leading attractor, we 
only take a middle point Vtm/^v — 0.3 as this turning point of the simulation. 

However when we estimate other theoretical values, it is impossible to adopt 
an exact manner. The quintessence has its energy density about p^s ~10^^^°Mp 
at the ending point 3 of the slow rolling, and the matter density is about pms — 
10~^^^Mp at this point. The evolution of the quintessence energy is like of the 
radiation, p^oc a~^, and in other hand pmoc for matter. When the cosmic scale 
factor expands to 04/03 = Pvs/Pms = 10^^, the quintessence energy will decrease 
faster and is approximated equal to the matter energy, the universe enters from 
the quintessence kinetic dominated phase into the matter leading attractor phase. 
The time is h/h ^ (04/03)^ = 10^^. Then it is easy to obtain pm4 — 10"^°"' 
and pc4 ^ 10-257]Vf4, The evolution of is still Eq.(|12D, one gets 04 = 99Mp and 

04~ All estimation is confirmed with the simulation well. 

The evolution of the matter leading attractor is like of matter, p^ = 3pm oc a~^, 
but the evolution of the curvature energy is slower, p^ oc a~^, thus when the cosmic 
scale factor expands 05/04 = Pm4/|Pc4| = 10^^, the curvature energy will dominated 
in the universe. The time at which the universe begins its curvature domination 
is ^5/^4 — (05/0.4)^/^ = 10^^. We choose fics — —1 as the definition of the turning 
point ts in the simulation. For same reason one gets 05 = 282Mp and 05~ lO^^^^M^ . 
Again all estimations are coincident with the simulation well. 

VIII. THE STOP OF THE EXPANSION AND THE START OF THE 

SHRINK OF THE UNIVERSE 

In the curvature energy dominated universe the curvature negative-energy will 
cancel strongly with the matter and quintessence energies, the Hubble parameter 
of the universe will become smaller and smaller. Ultimately, the Hubble parameter 
will become zero at the time point t^, and the universe stops its expansion, dg = 0, 
which is the balanced point. The time is not far from the time ^5, but very near, 
so that the values of the various parameters are not changed too much, og — 05, 
06 — 05, ^6 — ^5- We define the time difference t =t — t^. Through the observation 
of the various curves of neighbor of the balanced point of the simulation, we find 
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out the following expression is a good approximated solution for a narrow region 
|r| < O.Olte, 

= 06 + 06^ + ^06^^ a = aQ- iagr^. (14) 
Putting it into Eqs.(|^) and (^), we have 

= 2( \p,o\af - Mo^e-"'^^/^^'' - Pmoa^'), (15) 

0g = XM-^M^e-^'f'^'^^ (16) 

cie = Mp-=^(1 + 3A-2)-i(^|p,o|a6 ' - ^Pmo^ ')• (17) 

One can check easily that the theoretical estimation data obtained by using theses 
formulas are almost same with the simulation data. Let us look at the Hubbel 
constant, H = — a^^Ogr, we have H > when r < and H < when r > 0. The 
universe comes from the expansion to the shrinkage. It is very important that the 
sign of the Hubble parameter is changed, which turns the direction of the universe 
evolution. 



IX. THE BIG SHRINK OF THE QUINTESSENCE KINETIC ENERGY 

DOMINATION 

We have to note that at the balanced point, where Oe is at a maximum of the 
cosmic scale factor, the proportion of the quintessence kinetic energy is not small, 
Qk/^m = 0.94 (note that Qk,^m — oo at this time point). Once the universe 
comes into its contractive phase, the proportion of the quintessence kinetic energy 
will increase rapidly. A good approximation for this phase is to set = Pm = 
Pr = Pc = 0- This is a typical quintessence kinetic dominated model, the unique 
distinction is the sign of the Hubble parameter, which is now negative! But its 
evolution law is same, pk =0^/2 oc a~^, when the universe shrinks, the quintessence 
kinetic energy increases faster than any other energies. The just same behavior 
in the early expanding universe is used to construct a tracker solution ^^''l. The 
quintessence begins its speedy running at the very flat potential. Its running energy 
make universe to shrink more rapidly, and then the kinetic energy is higher. Using 
this simplification we obtain easily its solution for r > — from Eqs.(|^) and (^, 
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</. = </.6- V6Mp In (a/ae), a = a^^l - ^^^^M^^f/^ (18) 

If any unforeseen does not happen, the universe will shrink to zero size as = during 
the time interval rg = •y/2/30(^^Afp ~ 21^ from (|l^). In a view of the exponential 
time the universe shrinks to zero size almost suddenly. It is very different from the 
shrinkage of the matter dominated universe familiar for us. The time tg = i^^7 + ''"8 — 
lO^^^Mp"^ ~ 10^^^t„ is the life of a cycle of the universe in our model, which is an 
unimaginable very long time, but it is important that it is finite. 



X. THE SOME CONJUNCTION ABOUT A SMALL INITIAL 
CURVATURE ENERGY 

After rearrangement of the catenulate relations among the scale factors of the 
various evolution stages, we actually obtain a simple expression about the maximum 
cosmic scale factor for a cycle of the universe 

Pv2 

oe = -. ra2, (19) 

\Pc2\ 

i.e., it is mainly determined by the initial curvature energy. The formula is valid 
only if the ag is larger than 03, since only if the quintessence ends its slow rolling we 
are able to use this formula correctly. A focus question is why the curvature energy 
is so small in the initial condition? It involves a quite different stage of the universe 
evolution, i.e., the creation and inflation of the universe. We only can do a very 
simple discussions about this important problem since the major issue of this paper 
is about the present acceleration and the future constriction of the universe. 

Suppose that we have a inflation field $ which inflationary potential is the typical 
chaotic one'^1, i.e., V{^) = M2$72 (here the index "i" means "inflation"). Suppose 
that the inflaton begins its inflation at the initial field value $6 = 25Mp and ends 
its inflation at the final value $e = V^Mp with e($e) = 1 in our model. We see 
that this $6 is very far smaller than the critical value ^chaoUc > 300Mp at which 
the universe star its true chaotic inflation behavior'^^1 . The fluctuation 5^ ~ 10~^° 
of the CMBR observed by us will add a constraint on the important parameter of 
the inflationary model, i.e., the inflaton mass, then ~ lO^^^M^ is obtained from 
this constraint. Therefore the inflationary e-folding is = $^M~^/4 ^ 156, and 
we have = lO^^ab, where a^, is the size of the universe at its quantum creation. 
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In spite of this number is very large, it is actually very smaller than of the typical 
chaotic inflation. It is a good idea that the original universe comes from an four- 
dimensional instanton 5*^ with the radii in the Euclidean space (or to multiply 
directly another unvarying compact manifold in an extra dimension space to fit 
the need of the string (M-) theory^^^l). Only the closed universe {Kf. = 1) can be 
transformed from this instanton, and can supply a barrier to realize its quantum 
tunneling, which produces a model with the initial condition Hf, = and ${, = 0. 
This implies that the negative value of the curvature energy —p^ = K^a'"^ must 
equal to the infiaton field energy in Eq.(D at the born time of the universe^ ' 

Thus the we obtain ab — 3 x 10^M~^ and pcb = —^l ^ —10~^Mp, which is negative 
for the closed universe. At the ending of the inflation of the infiaton the curvature 
energy becomes pce — — 10~^^^Mp and the inflation field energy transform into the 
radiation energy p^e = P<i>e = = 10~^^Mp, this progress is the "reheating", 

which temperature is about Tg ~ pyj' ~ 10~^Mp. The BBN temperature is about 
To = lMeV~10~^^Mp. From the reheating to BBN, the universe expands by ao/oe = 
Te/To = 10^^, therefore we have the parameters at BBN time pco — — lO^^^^M^ and 
Pro — lO^^^Mp , which has been adopted by us in the section 3. Even if the reheating 
temperature is lower than Tg just given, we can consider a period of the matter 
dominated expansion of the coherent condensate state before the reheating and 
obtain a correct expanding ratio, here oq is only a naive estimation. The electro- weak 
interaction produces the dark matter, the baryogenesis produces the asymmetry of 
the baryonic matter, so that it is reasonable that the matter density is lower than 
the radiation one for many magnitude orders p^o — 10~^^Mp. If $5 is only small 
like 16Mp, the present scale factor has the same size of the magnitude order with 
the present event horizon of our universe. Therefore at least the initial value of the 
infiaton must be larger than 16Mp. Then we see that the initial value 25Mp of the 
infiaton set by us is not too higher than a necessary value for the realistic universe 
and far smaller than the chaotic critical value. Even so this has brought a very 
small curvature energy density, and the concrete value of pco taken here is only as 
an example. 

We would like say some words about the closed universe. If the universe was 
born by a quantum process, it is relatively easy to image that a finite size thing 



was produced from nothing (for example, instanton), but it is relatively difficult 
to image an infinite size thing was produced from nothing, in spite of this is a 
mysterious "quantum process" . Due to the curvature energy in the present universe 
is too small, Qcu = — 10~^^, it is very difficult to observe its existence. However 
we are not able to exclude this possibility of the closed universe, no matter how we 
improve on our observation precision. 

All just mentioned is high conjecture. We should distinguish the conjecture 
of this section with the reasoning in the previous sections. The curvature energy 
dependents on the inflaton initial value, which determine finally the life of a cycle 
of the universe. Anyhow, the cycle life is finite, but the number of the cycles may 
be infinite. 

XI. SOME CONJUNCTURE ABOUT THE BIG SHRINK 

We have owned the evolution rule of the quintessential kinetic energy shrink, 
i.e., (|18|) and pk oc a~^. It is convenience to define the ending time coordinate of 
the universe, i.e., x = tq — t > 0. Then we have a = Qq^^ 3 / 2^qM~^ xY^^ . As a 
good approximation, we take the balanced point 6 as the starting point, where the 
quintessence kinetic energy density is p^e — 10"^^'' M^, the matter density is pme — 
10~^^^Mp, the radiation density is prQ — lO^'^^^Mp and the curvature energy density 
is pcG — — lO^'^^^Mp. Since the quintessence has a large field value 06 — 282Mp 
in this time, and the quintessence potential energy is not able to increase as the 
universe shrinks, we neglect it hereafter. 

A question we have to discussed here is whether the matter has transformed into 
the radiation. This is possible. Since the proton will decay. However the baryonic 
matter is a small part of the whole matter, a large part of matter is the cold dark 
matter. The key problem is whether the cold dark matter will decay, many models 
think so. If a large part of the cold dark matter is the neutralino of the supersym- 
metric particles, and if the R parity is an absolute conserved quantum number, then 
the neutralino does not decay^^^'. To say the least, even if the neutralino decay, its 
final states must include the neutrinos, however maybe the neutrinos have a small 
non-zero mass. The the final state particles produced by the evaporation of the 
black holes are not the pure photons, they include the neutrinos at least. Therefore 
it is impossible that all matter can transform into the radiation. Of course even if 
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the matter transform into the radiation, it does not affect the quahtative conclusion 
about the big collapse of the universe in our model, only change some scenario, for 
example, the matter leading attr actor is replayed by the radiation leading attr actor. 

Now we would like to estimate a time point "gr" , i.e, "quantum gravity" , which 
will be explained soon. At this time point q the quintessence kinetic energy density 
arrive the Planck energy scale, i.e., pkq c± IM^. Note that this is very high density, 
you note that even if at the born moment of the universe and the begin of the 
inflation, the density is only 10~^Mp in the last section. It is easy to get ~ 
lO^^^ae ^ 10^°ao ^ 10^°a„, 0q ~ 290Mp and Xq ~ 1M~^ by using just mentioned 
formulae in this section. It is a very strange conclusion, since the energy density 
of the universe arrives the Planck energy scale at the Planck time until its ending, 
but the size of the universe at that time is larger than our present universe for 20 
magnitude orders! It is unimaginable. The quantum gravity has to play its role at 
this high energy scale in such wide space. How does this kind of universe with such 
high energy and such large space scope shrink further? 

We must note that at this time the quintessence mass is very small, but it is not 
zero, 

ml = V" ~ X^M-^M^ exp(-A0/Mp), for > 0, (21) 

and ~ lO^^^^M^ for 06- It is very possible that the rabid running of the 
quintessence transforms into its exiting states, i.e., particles due to its fluctuation. 
We do not know whether this transform is thermal. If its temperature is absolutely 
zero, its evolution is matter-like oc a~^. If its temperature is non-zero, its evo- 
lution may be radiation-like oc a""^ since its temperature will increase in the 
compression. We can re-estimate the time at which the density arrives the Planck 
scale, and get ag ~ 10~^°ao — 10~^°a„ for matter-like, and ~ ao — 10~^°a„ 
for radiation-like, both have Xg ~ 1M~^. Despite of this Planck scale universe 
is very smaller than the universe at BBN, but it is very larger than its born size 
Qb ~ 10^ M~^. According to this argument we think that the universe never compress 
to its original state at its birth. The new cycle must have its new property. 

A key problem is when, i.e., at what time, the quintessence kinetic energy trans- 
form to the particles? Maybe at the beginning of the universe big collapse, maybe 
at the beginning of the quintessence fast roUing. As mentioned, this does not affect 
our main conclusion, the matter dominated universe is almost same with the matter 
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leading attractor before the constriction. If the quintessence kinetic energy trans- 
forms to matter, the constriction of the universe is just hke the standard shrinkage 
of the mater closed universe. Same saying for the radiation-like. It is very possible 
that this transform happens at the contractive epoch, since the quintessence has a 
smaller mass and the universe at this epoch has very large deceleration which is far 
different from the inertia system. This epoch is a very different from other epochs 
with small deceleration. 

We see that the curvature does not play any key role when begin of shrinkage, 
its action is only to change the sign of the Hubble parameter. If there is a negative 
cosmological constant, the case is similar. But the quintessential kinetic energy 
exists still in the contractive stage even if for a negative cosmological constant, this 
is a new phenomena been worth to pay some attention on it. 

XII. CONCLUSION 

Our model realizes in a relative natural way an "universe" which is now acceler- 
ated expanding and will shrink ultimately. There is not the event horizon problem 
in our model. It seems that it escapes from the consequence of Ref. [6] . The way of 
shrinking is unexpected, i.e., the quintessence kinetic energy dominated, the velocity 
of the energy increasing is unexpectedly fast. Our model is semi- realistic since we 
obtain a time point "u" with almost same parameters of the present day universe, 
specially the universe age i„ and the acceleration parameter q^. Our model supply 
a cycle mechanism to realize the infinitely many universe based on the possible time 
cycle evolution, of course it is not excluded the possibihty of the existence of the 
other unconnected cychc closed universes. We do a computer simulation and the 
analytic analysis, both are well consistent. 

Our model is able to change the whole scenario about the origin and ruin, or 
recurrence of the cosmos. Our model begets some interesting problems to be studied 
worthily. 

An outstanding character of the "Niagara" potential is that the quintessence has 
a negative mass square in its slow rolling phase, i.e., its exited particle is a tackyon! 
We have 

m| = T/"~ -A2M-2Mo^exp(A(/)/Mp), for < -1, (22) 
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and ~ — lO^^^^M^ for cf) ~ —2. Whether this brings some observable effect is an 
very interesting problem. Kofman et. al. research the tackyon effect at the early 
inflation epoch of the universe'^^' . Maybe this tackyon effect can avoid a problem 
about the fifth force as an ordinary long range interactiont^^l . 

Other problem is the following. What is a more exact behavior about the evo- 
lution solution of our model. For example, how to transition exactly at the interim 
points 3 or 6? How to join exactly the approximation solutions Eqs.(^ and (|T^? 
How to join exactly the approximation solutions Eqs. (]T^) and (|T^? 

How to realize the rebound of the universe from a quintessential big collapse to 
a new big bang, or a new big infiation? 

What does the absolute high quintessence kinetic energy in the shrink stage 
transform into? Whether the quintessence kinetic energy can transform to matter 
or radiation in the contractive stage of the universe? Maybe the quintessence kinetic 
energy is also able to transform into a pseudo cosmological constant if we make a 
careful design for the quintessential potential, and then a radical different scenario, 
a part of defiation, will appear in the constriction of the universe. 

What is the origin of the "Niagara" potential? Whether we need to incorporate 
the early epoch tracker^-'^^] with it? 

Whether we still need the anthropic principle to adjust the theoretical parame- 
ters? We think so. 

If the universe will shrink, when it shrink to the size of the extra dimensions, 
what happens for the extra dimensions? 

I hope one can solve some of these interesting problems in the near future. 

We must supply a comprehensive understanding for universe from its birth to its 
termination. The unilateral quest of separate scenarios without an unitive viewpoint 
is not able for us to arrive a goal to uncover truly the riddle of our wonderful universe. 
The essence of the cosmological constant is various, itsits story is rich. 
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FIG. 1. The evolution of the acceleration parameter q by computer simulation. The 
first stage is the radiation dominated from to = 0.4 x 10^^ M^"^ to ti = 0.6 x 10^^ M^"^. 
The second stage is the matter dominated from ti to t2 = 0.86 x 10®'^ M^"^. The third 
stage is the quintessence potential energy dominated from t2 to = 3.2 x lO'^Uf-^ The 
fourth stage is the quintessence kinetic energy dominated from t^, to t4 = 1.6 x lO^^^Afp"^. 
The fifth stage is the matter leading attractor from ^4 to ^5 = 0.28 x lO-'^^^Afp"-'^. The 
sixth stage is the curvature energy dominated from to = 0.85 x lO^^^M^^^ The 
seventh stage is a big collapse with the quintessence kinetic energy dominated from tg to 
tg = 1.74 X W^'^M-^. Specially, tnow = tu = 1-56 x lO^OMp-^ = 13Gyr. 
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FIG. 2. The evolution of the cosmic scale factor a. Because its varying scope is too 
large, we have to use the double logarithm coordinates. In order to get a normalization, we 
add a minus one in it. You can see that the scale factor is shrinking in the stage 7. Some 
data are, ai = 1.2 x 10^ oq, a2 = 2.4 x 10^ ao, anow = = 4 x lO^ao, as = 1.7 x lO^^ao, 
04 = 1.8 X lO^'^ao, as = 3.1 x IO^^oq, oe = 4.5 x lO^^ao, aj = 0.27 x W^^ao, and ag = 0. 
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FIG. 3. The evolution of the quintessence field value (j). In order to look at a detail 
of the behavior in small field value, we use the logarithm coordinates. In order to avoid 
singularity for iniyial value — 2 , we have to add three in it. The curve is constructed 
by (p- of Eq.(8) for < and 0+ of Eq,(12) for > 0, there is a transition at time ts. 
Some data are: (^o = -2Mp, (/>! = -2Mp, 02 = -1.997Mp, (/)3 = 0.026Mp, 04 = 95Mp, 
05 = 280Mp, 06 = 281. 5Mp, and 07 = 289. IMp. 
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FIG. 4. The evolution of the velocity of the quintessence, d(j)/dt. This figure includes 
more information. You can magnify it in your computer to watch a detail for the main 
curve, specially, points 3 and 5 — 8. Some data are: 0o = 0, 0i = 0.63 x 10^^ ''M^, 



^2 = 0.8 X IQ-^'^M^, (j)3 = 0.32 X IQ-^^M^, 04 = 0.63 x IQ-^^^M^, 



0.4 X IQ-^^^M^, 



0.21 X IQ-^^^M^, and (pj 



10 ^''^Mp. It is notable that the point 6 is not at the 



minimum of the curve, but is at its right side, therefore, d?(j)/dt'^ > 0, see Eq.(16). 



